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ABSTRACT 

A recent analysis of the Einstein sample of early-type galaxies has revealed 
that at any fixed optical luminosity Lb SO galaxies have lower mean X-ray 
luminosity L x per unit L B than ellipticals. Following a previous analytical 
investigation of this problem (Ciotti & Pellegrini 1996), we have performed 2D 
numerical simulations of the gas flows inside SO galaxies in order to ascertain the 
effectiveness of rotation and/or galaxy flattening in reducing the Lx/Lb ratio. 

The flow in models without Supernova (SNIa) heating is considerably 
ordered, and essentially all the gas lost by the stars is cooled and accumulated 
in the galaxy center. If rotation is present, the cold material settles in a disk 
on the galactic equatorial plane. Models with a time decreasing SNIa heating 
host gas flows that can be much more complex. After an initial wind phase, 
gas flows in energetically strongly bound galaxies tend to reverse to inflows. 
This occurs in the polar regions, while the disk is still in the outflow phase. 
In this phase of strong decoupling, cold filaments are created at the interface 
between inflowing and outflowing gas. Models with more realistic values of the 
dynamical quantities are preferentially found in the wind phase with respect 
to their spherical counterparts of equal L B . The resulting L x of this class of 
models is lower than in spherical models with the same L B and SNIa heating. 
At variance with cooling flow models, rotation is shown to have only a marginal 
effect in this reduction, while the flattening is one of the driving parameters for 
such underluminosity, in accordance with the analytical investigation. 

Subject headings: Galaxies: elliptical and lenticular, cD - Galaxies: ISM - 
Galaxies: structure - X-rays: galaxies 
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1. Introduction 

It is well known that normal early-type galaxies are X-ray emitters, with 0.2-4 keV 
luminosities ranging from ~ 10 40 to ~ 10 43 erg s _1 . The X-ray luminosity Lx is found to 
correlate with the blue luminosity Lb (Lx oc L b ' 0±0 ' 2 ), but there is a large scatter of roughly 
two orders of magnitude in Lx at any fixed Lb > 3 x 1O 1O L flFabbiano 1989| ; Fabbiano, 
Kim, & Trinchieri 1992). 

The trend and the scatter in the Lx — Lb diagram have been successfully reproduced 
by the time-dependent, ID hydrodynamical models of Ciotti et al. (1991, hereafter CDPR). 
In these models the gas lost by the evolving stars is heated to X-ray temperatures by stellar 
random motions and SNIa explosions. Given the assumed time dependence of both mass 
loss rate (oc t~ L36 ) and SNIa rate (oc t -1 - 5 ), the CDPR models can evolve through up to 
three different dynamical phases, i.e., from winds to subsonic outflows to inflows (in the 
following we call this kind of evolution the WOI sequence). CDPR introduced the global 
parameter x m order to predict the dynamical status of the gas flows in non-rotating, 
spherical galaxies. This parameter was defined as the ratio between the power required to 
extract steadily the gas supplied by the mass loss of the evolving stellar population from 
the galaxy potential well, L g , and the power supplied to the gas by the thermalization of 
the stellar random motions and by the supernova heating (L a and Lsn, respectively). Due 
to the temporal evolution of the various source terms, % is an increasing function of time. 
Numerical simulations showed the correspondence between x < 1 an d the presence of a 
wind phase (i.e., of a low Lx), and the onset of a high Lx soon after x becomes > 1. 

A recent analysis of the complete Einstein sample of early-type galaxies showed that 
on average SO galaxies have lower Lx and Lx/Lb than ellipticals (Es) (Eskridge, Fabbiano, 
& Kim 1995a,b). A possible explanation of this fact could be that SOs and non-spherical 
Es are less able to retain hot gaseous haloes than rounder systems of the same Lb- This 
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could be due to a lower gravitational energy, because of their higher rotation rate which 
decreases the effective potential, or their different mass distribution. A theoretical analysis 
of this problem - based on the global energetics of the X-ray emitting gas - was carried out 
by Ciotti & Pellegrini (1996, hereafter CP). They generalized the global parameter x f° r 
flows in flat and rotating two-component galaxy models. In particular, they added to x the 
contribution of the ordered rotation of the galaxy stellar component, L rot , showing that, at 
least from a global point of view, rotation alone cannot change an energetically bound gas 
flow to an unbound one (see §2.3). On the contrary, using two-component Miyamoto-Nagai 
models (1975, hereafter MN), and under the assumption of an identical rate of SNIa per 
unit L B in SOs and Es, CP showed that flat galaxies can have a much lower x than rounder 
ones of the same luminosity, due essentially to their flattening, even in the presence of 
massive spherical dark haloes. More precisely, they showed how a realistic flattening can 
produce a decrease of L g of ~ 20%, a result which can be sufficient to change an inflow 
into a wind. So, the effect of the flattening was suggested as a possible explanation for the 
SO underluminosity. Incidentally, note that the decrease in the gravitational energy that 
is released by a stationary inflow at the bottom of the galactic potential well, associated 
to the decrease of L g , does not lower significantly the Lx of pure cooling flow models (cfr. 
Fig.l in CDPR). 

Despite the possible agreement between the observations and the results of CP, this 
purely energetical approach can be seriously questioned considering the complexity of 2D 
gas flows and the fact that x is a global quantity. In fact decoupled flows - in which at the 
same time different parts of the galaxy are in inflow and outflow - may in principle develop 
in the CDPR scenario in such a way that Lx remains high even if % < 1; if this is the case, 
X < 1 cannot be used as a reliable indicator for a low X-ray luminosity, and the analysis of 
CP is seriously questioned. 
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2D numerical simulations of gas flows in early-type galaxies were actually carried 
out by Kley & Mathews (1995, hereafter KM) and Brighenti & Mathews (1996, hereafter 
BM). These simulations, however, considered only pure cooling flows for which the flow 
decoupling is not expected. KM considered rotating spherical King models plus a spherical 
quasi-isothermal dark halo, without SNIa heating. Successively BM considered oblate 
galaxies with different axial ratios, sustained by various amounts of ordered and disordered 
kinetic energies. They assumed a SNIa rate oc t _1 , with a present day value of one-tenth of 
that proposed by Tammann (1982). Their models show that the gas flow and the associated 
Lx are only slightly dependent on the flattening of the galaxy (in agreement with the CP 
analysis applied to the inflow case), and that L x is reduced by increasing the amount of 
ordered rotation. These 2D hydrodynamical models do not solve problems such as the 
disposal of excessive amounts of cold gas (a problem already encountered in ID pure cooling 
flow models) on an equatorial disk, which is actually unobserved, and the excessively flat 
shape of the X-ray isophotes when compared with the observations. 

Motivated by the given arguments, we adopt the same input physics used for the 
hydrodynamical equations in CDPR, in order to explore whether the WOI sequence is still 
achieved in flat galaxies, how such a sequence is influenced by galaxy flattening and rotation, 
how important the decoupling in the flow is, which observational consequences are predicted 
in this case, and finally whether the energetical explanation of CP is exhaustive. This is 
accomplished with two different classes of models. In the first we compare the behavior of 
the gas flow in spherical and flat galaxies that are energetically globally equivalent (both 
with and without SNIa heating). In the second class of simulations we investigate the 
gas flows in realistic galaxy models of SOs, also in order to explore the possibility of flow 
decoupling in galaxies with x < 1 an d its effects on Lx- 

We thus study the evolution of gas flows in the galaxy models constructed and described 
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in detail by CP, i.e., axisymmetric stellar and dark mass distributions with the MN shape, 
and in which the internal dynamics is determined by the solution of the associated Jeans 
equations, allowing for variable amounts of ordered rotation and velocity dispersion. 

In Section 2 we briefly describe the models used in the simulations. In Section 3 
the results obtained for models with different amounts of ordered rotation, SN heating, 
flattening, and characteristic scale-length are presented. A discussion, considering also some 
recent observational results is given in Section 4, and in Section 5 the main conclusions are 
summarized. 



2. The models 

2.1. The galaxy density distribution 

The two-component mass models used in our simulations and their internal dynamics 
are fully discussed in CP, which includes the explicit analytical formulae for the various 
terms in the virial theorem and for \ are given. Here only the basic formulae are reported. 
The stellar density distribution is given by 

(M.*b 2 \ aR 2 + (a + 3() (a + C) 2 m 



and its potential is 



4tt J ( 3 [R 2 + (a + C) 2 ] 5/2 ' 



= == ■ (2) 

/ R 2 + (a + () 2 



where ( = V z 2 + 6 2 , and (R, ip, z) are the cylindrical coordinates (MN)[]. The dark matter 
halo is described by another MN density distribution characterized by Mh, ah and The 



1 When a = the density distribution (1) reduces to the Plummer (1911) model, and b 
becomes its effective radius. So, b can be roughly considered a characteristic scale-length, 
while the galaxy flattening is determined by the dimensionless ratio s = a/b. 
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associated Jeans equations are solved in the total potential $ tot = + $ h , and under 
the assumption that the distribution function depends only on the specific energy and the 
angular momentum component along the z-axis (see, e.g., Binney & Tremaine 1987, p. 197). 
With this assumption the radial and axial velocity dispersions are equal, i.e., <7r = er z = a, 
and the only non-zero streaming motion can be in the azimuthal direction (<p) . The intrinsic 
degeneracy in the azimuthal velocity component is removed adopting the Satoh (1980) 
/c-decomposition, i.e., 

v? = k 2 (^-o\ (3) 

and 

al = ^-^ = k 2 a 2 + (l-k 2 )^ (4) 

with < k < 1, where is the modulus of the ordered rotational velocity of the galactic 
stellar component. In the case k = no net rotation is present and all the flattening is due 
to a 2 ;k = l corresponds to the isotropic rotator. The analytical solutions of the Jeans 
equations can be obtained for 6 h = b. In simple case of a h = a the formulae for the velocity 
components are expressed in terms of elementary functions and are given in CP together 
with their edge-on projection. 



2.2. The hydrodynamical equations 

The equations describing the gas flows in our simulations are 

dp 

— + V • (pv) = ap„ (5) 

dm _ 

— + V-(m®v) = -( 7 -l)V£ + gp + ap,v, (6) 
at 

_ + V • (Ev) = -( 7 - 1) EV ■ v - L + ap*(e + - ||v - v*|| 2 ), (7) 

where v is the fluid velocity, and p, m, and E are the mass, the momentum, and the 
internal energy of the gas per unit volume, respectively. The source terms on the r.h.s. of 
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equations (5)-(7) describe the injection of mass, momentum and energy in the gas due to 
the mass return and energy input from the stars. a(t) = a*(t) + q;sn(*) is the specific mass 
return rate from stars and SNIa respectively, with a* oc t~ im and «sn oc t~ L5 ; e is the 
injection energy per unit mass due to the stellar random motions and to SNIa explosions. 

is the bulk velocity of the stars, whose modulus is given by equation (3), g = — V$ to t 
is the gravitational acceleration due to the total mass distribution, and 7 = 5/3. Finally, 
L = n c n p A(T) is the cooling rate per unit volume, where for the cooling law A(T) we follow 
the prescription of Mathews & Bregman (1978). For more details on the exact form of 
«*(£), o:sn(0 an d e o, see CDPR and KM. 

To integrate the set of equations we used a second-order, upwind, 2D extension of the 
numerical code described in CDPR, coupled with a staggered, spherical, Eulerian grid. 
Given the symmetry of the problem, the grid covers polar angles 0° < 9 < 90° with 40 
equally-spaced angular zones, and reflecting boundary conditions are applied at 9 = 0° and 
9 = 90°. The radial coordinate covers the range < r = y/R 2 + z 2 < 100 kpc with 80 
logarithmically spaced zones, and a free outflow from the grid is allowed at the outer edge. 
At the inner edge reflecting boundary conditions are set. We assume the model galaxies to 
be initially devoid of gas due to the previous activity of the Type II SNe. A discussion on 
the reliability and implications of this assumption is given in CDPR. 



The main ingredients of the global parameter x are presented here. The SNIa heating 
evolution is parameterized as L SN = 10 51 i?sN-^B erg s _1 , with 



2.3. The x parameter 
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where h = ^/lOOkms" 1 Mpc" 1 , and 1 SNU = 1 SN per century per 1O 1O L (cfr. CDPR). 
When t = 15 Gyr, L B = 10 10 L o and $sn = 1, the standard SNIa rate is recovered 
(Tammann 1982)0. 

The power L g required to extract steadily from the galaxy the mass losses from the 
aging stellar population is 

L s = a(t) J p*|$tot|dV. (9) 
The thermalization of the stellar random motions heats the gas at a power L a given by 

r _ «(*) 



2 



J p*(2a 2 + ffJ)dV, (10) 



while the power L rot related to the ordered fraction of the kinetic energy is given by 



J rot 



J pjtfdV. (11) 



Unfortunately this last term cannot be inserted in the \ definition in a simple way. 
In fact, two opposite and extreme scenarios can be assumed: one in which also L rot is 
thermalized (contributing to the gas heating together with L a and Lsn) ; and another in 
which the only effect of L rot is to reduce the depth of the effective galactic potential well. 
As a consequence, the global parameter \ is defined as 

x= T /;:]f rot u , (12) 

L SN + L a + (1 - 7)L rot 
with < 7 < 1. Note that if 7 = 0, i.e., in the case of complete thermalization, the virial 
theorem implies that L a + L rot depends only on the galaxy structure and % is independent 
of the amount of ordered rotation (in our case of the value of the parameter k). When 
instead 7 = 1, i.e., in the case of cold rotation, the effect of rotation is maximum; CP 



2 Note that in our models the assumed Lb is independent of h. As a consequence Lsn 



depends on h. In the paper we adopt h = 1/2 
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showed that also in this case, for realistic galaxy models the differences in the value of x, 
for a null and a maximum ordered velocity configuration, are small, usually less than 10%. 
Moreover whenever 7 > the introduction of rotation decreases a x < 1 and increases a 
X > 1. In conclusion, from the CP analysis it results that rotation alone cannot invert the 
global energy balance of a gas flow, but can help the occurrence of the wind when x < 1- 
Due to the small dependence of x 011 7, through the paper we assume 7 = 0. 

3. The results 

3.1. Model Al: $ S n = 0.3 and k — 1 

We describe here in detail the main characteristics of the flow hosted by a galaxy model 
in which ordered rotation (k = 1) and SNIa heating in accordance with the latest estimates 
of Cappellaro et al. (1997, $ S n = 0.3) are present. We assume M* = 2.75 x 10 11 M , 
M^/Lb = 5.5, so that Lb = 5 x 10 10 L Q , and a = b = 3 kpc. The dark halo has the same 
profile of the stellar distribution, i.e., = a and &h = b, and Mh/M* = 2 (see Table 1). 
The three-dimensional central velocity dispersion is a c = 330 km s -1 , and the maximum of 
the ordered rotational velocity is v Iot = 360 km s -1 . This set of parameters is such that 
X = 1 at the present epoch, in order to compare as much as possible the present results 
with those obtained for the so-called reference model discussed in CDPR (hereafter KRM, 
Xkrm = 1-1). L B is also the same of the KRM, so that the mass return rate of the two 
models is exactly the same. Note how it turns out that, in order to keep x — 1> model 
Al must have very high stellar a c and t> rot . This is not a surprise: in fact, as already 
pointed out by CP, in a realistic flat model with a non-negligible $sn the present-epoch x 
is substantially lower than unity. The same comment applies to models A1-A4. 
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3.1.1. Hydrodynamics 

As a consequence of the early high SNIa energy injection, due to the assumed $sn 
value and «sn oc t _L5 , at the very beginning the flow is in the wind supersonic phase, i.e., 
the radial velocity u r (R,z) > v somid (R,z) throughout the galaxy (Fig.l, top left panel). 

As time increases and the energetic input decreases according to equation (8), the 
sonic surface moves outward, with u r decreasing earlier in the central regions of the galaxy, 
and the gas density increases. Finally, due to the rapidly increasing radiative losses (the 
so called cooling-catastrophe), at t cc ~ 2 Gyr the fluid starts to revert its motion in a 
conical region containing the polar axis. We call this phenomenon flow inversion, and the 
subsequent flow decoupled inflow. The most striking feature of this evolutionary phase is 
the formation of cold, dense filaments at the boundary where two streams of inflowing and 
outflowing gas converge, as can be seen in Fig.l (top right panel). This filamentary phase, 
which in model Al suddenly develops at t m ~ t cc , is rather short (^2 Gyr). Later on, no 
more filaments are generated, and the previous ones are accreted at the galactic center. 
However, regions of gas colder than X-ray temperature and with a low density contrast are 
present for a much longer time. During this evolutionary phase two qualitatively different 
regimes are present: an inner zone (r<15 kpc) in which the fluid is rather turbulent, and an 
outer zone where the velocity field is smooth (Fig.l, bottom left panel). The inflow region 
becomes progressively larger and larger. At t — 15 Gyr the gas is still outflowing in the 
outskirts of the disk. The flow in the rest of the galaxy is in a turbulent state, characterized 
by very subsonic velocities (||v|| < 5 10 6 cms _1 , Fig.l, bottom right panel). 

In order to understand better the structure of the cold filaments, we performed a 
high-resolution simulation (with 300 radial meshes and 140 angular meshes covering the 
same physical space as the normal resolution simulations). Because of the severe limitation 
due to the azimuthal time-step, we could not extend the simulation over 15 Gyr. We thus 
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remapped onto the finer grid the hydrodynamical variables computed on the coarse grid at 
a time immediately before the filament formation, and we carried on the simulation for the 
entire lifetime of the filament. Fig. 2 shows the temperature and velocity field at t — 2 Gyr 
Density and temperature vary along the filament. For r<20 kpc, p&\ ~ 6.3 10~ 27 g cm -3 , 
and the density contrast between the cold filaments and the surrounding density is of the 
order of 30. The mean temperature is Tgi ~ 10 4 K (the lowest allowed temperature in our 
simulations). At large r the density decreases and the temperature increases, but never 
reaches X-ray temperatures. We did not made convergence tests, so it is plausible that the 
filaments are denser and thinner. We stress that in our 2D simulations these filaments are 
in fact surfaces. This fact should be kept in mind in all the following discussions^} 

The cold gas accreted during the entire lifetime of the galaxy forms a small disk, 
whose radius oscillates around a mean value of 3 kpc. The total mass of the cold disk is 
1.5 1O 1O M , while the hot X-ray gas is 6.3 10 8 M Q (Fig.3b, dotted line 1), much less than 
the hot gas in the equivalent KRM in CDPR (~ 5.4 10 9 M ). 

3.1.2. Energetics 

The comparison with the KRM is instructive. Also for this spherically symmetric 
model the present day x — 1> but t cc ~ 10 Gyr, much later than in model Al. From Fig. 3d 
we note that the time evolution of Lx (dotted line 1) is qualitatively similar to that of 
spherical models. After an initial decrease during the wind phase, Lx increases up to a 
maximum value at t — 2.5 Gyr. Later on, Lx decreases again, in pace with the decreasing 

3 The formation of such filaments seems to be a rather general feature in decoupled 
flows such as this. Melia, Zylstra, & Fryxell (1991), for instance, find astonishingly similar 
structures in their studies on accretion disk coronae around black holes. 
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mass return from the stars. When L x reaches its maximum, a non-negligible fraction of gas 
is still outflowing from the galaxy, as in the ID models. Note that at the flow inversion the 
energy budget (x < 1) indicates that the SNIa energy injection would be able to sustain a 
global outflow. This means that, because of the strong decoupling of the flow, in this model 
the global energetic argument (i.e., the x value), is not a useful indicator of the dynamical 
state of the gas. 

A second important difference between spherical models and model Al must be 
stressed: in the latter, Lx is a factor of four fainter than Lsn (Fig-3d, solid line), at variance 
with ID models in which Lx may be even greater than Lsn in the inflow phase. The 
fainter Lx of model Al is due to its lower content of hot gas with respect to the KRM (see 
§3.1.1) . This lower fraction of hot gas is due to the higher efficiency of the cooling near 
the equatorial plane of model Al, a natural consequence of a higher gas density in this 
region, which in turn is due to the higher local stellar density (the effect of rotation will be 
discussed in §3.1.4) 

We note however that, although the hot gas mass in model Al is lower by about a 
factor ten with respect to the KRM, L x does not diminish by the factor one hundred that 
one could expect as a rule of thumb (the cooling per unit volume being oc p 2 ): the much 
lower than expected reduction of Lx is thus due to different hot gas distributions between 
the two models. 

We finally point out that the noisy temporal behavior of L x is a consequence of the 
density and temperature fluctuations associated to large eddies in the flow. 
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3.1.3. Surface Brightness 

The X-ray surface brightness (E x ) evolution of model Al is shown in Fig.4 at the 
same epochs of the shots in Fig.l. The optical isophotes (£*) are also superimposed for 
comparison, and the galaxy is seen perfectly edge-on. 

During the wind phase, the X-ray isophotes are peanut- shaped, with regular lobes 
aligned with the galactic polar direction (Fig. 4a). The elongation of the isophotes decreases 
steadily with time, they become rounder and rounder, and, in the inner regions, eventually 
flatter than the optical isophotes (Fig. 4b, c). The isophotes of E x at the present epoch are 
rather flat in the central regions, and become progressively circular moving outward. This 
is due to the well known fact that the emitting gas accommodates in a potential well that 
is more spherical at large r. The irregularities in the isophotal shapes are clear evidence of 
the inner turbulence (Fig.4d). 

Note that once the cold filaments appear (Fig.4b,c), E x becomes very disturbed: we 
stress however that such an effect is magnified in our simulations by the imposed geometry, 
in which the filaments are in fact surfaces. We also computed E x under different viewing 
angles. It turns out that at early epochs the X-ray isophotes appear still elongated along 
the optical minor axis, even for values of the viewing angle as low as 45° (90°=edge-on). 
Instead, at late epochs, E* and E x become very similar as soon as the inclination angle is 
slightly lower than 90°. 

3.1.4. Model A2: $ S n = 0.3 and k = 

We describe here a model very similar to Al, but with k — 0, all the galaxy flattening 
being due to (see Table 1). For this model the value of x is independent of any 
assumption on the thermalization of ordered motions [see eq. (12) with 7 = 0], and is equal 
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to unity. 

From Fig. 3b it is apparent how this non- rotating model has a higher total gas mass 
(dashed line 2) than the rotating model Al (dashed line 1). This is because rotation - if 
not entirelly thermalized - decreases a x < 1 ( see §2-3), and in fact the largest difference 
between dashed line 1 and dashed line 2 takes place during the initial wind phase, when 
X is substantially lower than unity. Also the hot gas mass of model A2 (dotted line 2) is 
higher than that of model Al (dotted line 1): this is due to the combined effect of rotation, 
already discussed, and of the presence of a slightly more massive cold disk in model Al, a 
consequence of angular momentum conservation. This effect is quite small in the models 
heated by SNIa, but much more important in the cooling-flow models (a more detailed 
discussion is postponed to §§3.2.2-3.2.3). 

The main features of the temporal evolution of L x are similar for models Al and A2 
(Fig. 3d dotted line 2). We note however that the L x of model A2 is slightly higher than 
that of model Al, in accordance with its higher content of hot gas. The flow dynamical 
evolution of model A2 is also quite similar to that of model Al. The transient phase 
of filamentary instabilities - although to a minor extent - is still present, and develops 
approximately at the same time as in model Al. The absence of rotation does not affect the 
presence of oscillations in the temporal evolution of Lx, thus indicating that the turbulence 
is mainly due to a decoupled inversion of the gas flow, rather than to rotation. Finally, S x 
evolves mainly as in model Al, although the X-ray isophotes appear more circular in this 
non-rotating case. 
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3.1.5. Models A3 and A^: $ SN = 0.3 and k = I, but different a/b's 



We focus here on the relationship between the galaxy's flattening and the development 
of the cold instabilities. To investigate this point we run two models similar to model Al, 
except for the ratio a/b (see Table 1). A quantitative measure of the flattening of a MN 
model in its central regions is obtained with a series expansion of eq. (1). Retaining the 
quadratic terms the resulting isodensity surfaces are similar ellipsoids whose ellipticity is 



and s = a/b. Thus s = 1 adopted in models Al and A2 corresponds to S p = 0.45. As 
a consequence, with s = 0.7 (and £ p = 0.36) model A3 is slightly more spherical than 
model Al, and its x — 1-2 is higher, as expected. The result is that the flow inversion 
occurs earlier, and the filamentary phase is already completely developed at = 1.4 Gyr. 
On the contrary model A4, with s = 1.3 and and £ p = 0.52 is flatter. Its x — 0-94, and 
the filamentary phase starts now at tm = 3 Gyr. For a short time this model shows the 
contemporary presence of two cold structures similar to that in Fig. 2. These two models 
show how small changes in the galaxy flattening can change the flow evolution in the CDPR 



A very different flow evolution is obtained by suppressing the SNIa heating in a model 
equal to Al: now x — 4 since the beginning, and the behaviour is that of a pure cooling 
flow. The velocity field is very smooth, thus showing that SNIa's play a fundamental role 
in maintaining a boiling hot gaseous halo. In this case we do not obtain the transient 
filamentary structure, and Lx decreases smoothly. The present day Lx is more than one 
order of magnitude lower than in the previous SNIa heated cases (Fig. 3c, dashed line), 




(13) 



scenario. 



3.2. Model Bl: $ S n = and k = 1 
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because there is considerably less hot gas in the galaxy (Fig.3a, dashed line). The cooling 
gas accumulates now on an extended cold disk of ~ 10 11 M , approximately ten times more 
massive than that of model Al. Such a disk, present also in the simulations of BM and KM, 
is not actually observed in real galaxies. The X-ray isophotes of model Bl closely follow the 
optical ones. 

3.2.1. Model B2: tf SN = and k = 

Model B2 is analogous to model Bl, but no ordered rotation is applied. In this case Lx 
is higher than in model Bl at every time (Fig.3c, dotted line), just due to the larger amount 
of hot gas (Fig.3a, dotted line). At first sight it could seem contradictory that the model 
with the lower L x (Bl) has cooled a larger amount of gas, i.e., one could expect a higher 
Lx, at least in the past. However, although more mass cools per unit time in the rotating 
model, the lack of thermalization of the gravitational energy of the gas which stops on the 
rotating cold disk instead of falling into the galaxy center produces a net quenching in the 
global Lx (as pointed out by BM). The Ex isophotes in the non-rotating case are more 
concentrated towards the center and rounder than in the rotating found by BM. 

It is of interest to compare our cooling flow models with those of BM. These authors 
find a dramatic reduction of L x when introducing rotation, but at variance with us, the hot 
gas mass in their rotating models decreases only slightly compared to the non-rotating ones. 
Thus, the lower luminosity in their rotating models is due essentially to a redistribution 
of the hot gas mass, rather than to its scarcity as in ours. The origin of this different 
behaviour is in the higher flattening of the models used here (i.e., in the presence of the disk 
in the MN density distribution). In fact, in our models a larger fraction of gas is produced 
at \ow-z over the galactic equatorial plane. When the model is rotating this gas falls on the 
disk approximately at the same distance from the center where it was formed, bacause of 
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angular momentum conservation. So, on one hand it releases a much lower gravitational 
energy, and, on the other, the hot gas distribution is changed only slightly with respect to 
the non-rotating case. The accumulation of the gas on the disk also significantly increases 
its cooling, and the final result is a reduction of the hot gas mass. Note that also in the BM 
models the amount of cold gas increases with an increasing ellipticity of the parent galaxy. 

3.3. Models CI and Dl: $ SN = 0.3 and k = 1 

As mentioneded in the Introduction, we also ran models which are similar in many 
aspects to Al, but with larger 6's, in order to have realistic values for the velocity dispersion 
(see Table 1). Since p*(0) = M*(3 + s)/47r6 3 (l + s) 3 , increasing b at fixed s produces a 
reduction in the gas density and thus in the cooling, because of the lower stellar density. 
At the same time, because $*(0) = GM*/b{\ + s), L g is reduced too, favoring the escape of 
the gas. 

Model CI (6 = 4 kpc, er c = 290 km/s, i> ro t = 309 km/s, x — 0.85) develops a decoupled 
inflow at 5.5 Gyr, as can be seen from the evolution of Lx (Fig.3d, heavy dotted line 3). 
The most striking feature is again the formation of cold filaments (Fig.5). Surprisingly 
enough, these filaments do not form at the moment of flow inversion (t cc ~ 5.5 Gyr), but 
quite later (t m ~ 10.5 Gyr). At variance with model Al, they form close to the galactic 
equatorial plane (cfr. Fig. 2). Once on the plane, they slowly drift towards the center and 
accrete on a small cold disk (r ~ 4 kpc, M disk ~ 5 10 9 M ). This dynamical phase is 
characterized by a conspicuous and steady decrease of Lx, a consequence of the decrease 
of the hot gas mass (Fig.3b, heavy dotted line 3). While this effect is certainly real, it is 
magnified in our simulations for two reasons. First, as already pointed out, the filaments 
are in fact cylindrical surfaces; second, as usual in 2D simulations, reflecting boundary 
conditions on the disk introduce non-physical effects favoring cold material accumulation. 
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The S x evolution is similar to that of model Al, but postponed in time. 

For model Dl (6 = 5 kpc, a c = 258 km/s, v TOt = 277 km/s, x — 0-7), the resulting flow 
is a pure wind all over the galaxy, lasting for all the Hubble time without any decoupling. 
All the gas is hot, and its mass - as is common in the wind phase - decreases steadily and 
is very low (Fig. 3b, heavy line 4). L x is also very low and decreases proportionally to the 
gas mass (Fig. 3d, heavy dotted line 4). S x is similar to that shown in Fig. 4a. This run 
clearly shows how in a realistic flat model, even with a low i?sn, the flow can be in the wind 
phase during its entire evolution. 

3.3.1. Models C2 and D2: $ S n = 0.3 and k = I, but different a/b's 

We now briefly describe the behaviour of two models similar to CI and Dl, but with 
different flattening (see Table 1). The first (C2, s = 1.3, S p = 0.52, x = 0.75) is flatter than 
model CI, and consequently has a lower x- The resulting flow is a wind during all the time. 
This model (with 6 = 4 kpc) is energetically very similar to Dl (s = 1 and 6 = 5 kpc), as 
testified by their x values, and its flow evolution is also similar. 

Model D2 (6 = 5 kpc, s = 0.6, S p = 0.33, x = 0.86) is obtained from Dl by decreasing 
its flattening. As a consequence the gas is more bound as testified by a x similar to that 
of model CI (s = 1 and 6 = 4 kpc). The resulting evolution is similar to that of model 
CI, with the flow inversion and the following development of a decoupled inflow around 
t cc ~ 5.5 Gyr. However, in this case no cold filaments are developed. In conclusion, a very 
small change in the flattening can compensate for a variation of the characteristic scale 
length, showing again the strong sensitivity of models to variations of galactic structural 
parameters in the CDPR scenario. 
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4. Discussion 

Although the number of models we ran is not sufficient for an accurate statistical and 
quantitative analysis of the gas flow properties in flat and rotating galaxies, a qualitative 
comparison with the spherical models discussed in CDPR and with cooling flow models 
is possible; it is also interesting to compare our results with recent observational findings 
(Pellegrini, Held & Ciotti 1997). 

In order to compare the behavior of SO models with the equienergetical spherical 
models discussed in CDPR, we performed calculations for model galaxies with % = 1 at the 
present epoch (models A1-A2). We stress again that while x = 1 is a quite natural product 
of realistic spherical models, it turns out to be rather extreme for flat models. In any case, 
in flat models a clear X-ray underluminosity with respect to their spherical counterpart 
is found, mainly due to a drastic reduction of the hot gas mass. The amount of galactic 
rotation contributes only slightly to such a reduction. 

We also ran a few models with -# S n = (models B1-B2). These have a remarkably low 
L x , in contrast with spherical cooling flows which in general show excessively high X-ray 
luminosities. Although the model without rotation displays a present day L x similar to 
that of the rotating model, it was much more luminous in the past. As a consequence, with 
different choices of the structural and dynamical parameters, a large difference in the Lx 
can be obtained. Such a possibility is supported also by the models of BM, where different 
degrees of rotation lead to a spread in L x . 

When the galaxy structure and dynamics are more similar to those of observed galaxies, 
and SNIa's are present, x is l° w > an d the galaxy is permanently in a global wind phase or 
loses a large amount of gas during a Hubble time; the flow decoupling takes place at late 
times, even for -#sn as low as 0.3. The qualitative predictions made in CP are thus correct: 
it is much easier to extract gas from a flat, realistic galaxy model, than from a rounder 
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system of the same L B . L x may span all the values in the range between that of a pure 
wind and approximately Lsn by, for instance, varying b or the galaxy flattening, similarly to 
what happens for spherical models in the CDPR scenario for small variations of the model 
parameters. 

Moreover, the X-ray underluminosity of flat galaxies can naturally be accounted for 
because their Lx never reaches Lsn, even in the inflow phase. A possible explanation of 
the observations could be that at low L B , where the CDPR scenario predicts a significant 
number of objects to have an energy budget sufficient to unbind the gas, flat galaxies are 
preferentially found in the wind phase with respect to spherical galaxies of the same optical 
luminosity, while at high Lb, where a significant number of galaxies contain inflows, the 
underluminosity of flat galaxies could be due to the underluminosity of their inflow phase, 
as found (albeit for different reasons) in our models and by BM. 

We computed also a few models (not shown here) with a different flattening of the dark 
matter distribution with respect to the stellar component. The results are qualitatively 
similar to those discussed above. As a general rule, for a given stellar distribution, a rounder 
halo increases x an d reduces tm- For example, a model similar to Al but with a spherical 
halo (a/6)h = 0, has x = 114 and tm ~ t cc = 1 Gyr; a model similar to C2 but with a 
rounder halo {a/b\ = 0.6, has x — 0-8 an d im — 6.2 Gyr. As expected, as the flattening of 
the dark halo increases, x decreases and tm increases: for example, a model similar to D2 
but with (a) = 1.3 always remains in the wind phase, in accordance with the lower value 
of x = 0.79. 

We conclude that, on the basis of L x alone, it is not possible to decide whether the 
CDPR scenario or the cooling flow scenario is more suitable to describe the flows in SOs 
and Es: both can reproduce the X-ray underluminosity of flat galaxies. The scatter in 
the Lx — L-q diagram (at a fixed optical luminosity) is mainly produced by the galactic 



-22 - 



flattening (and rotation to a minor extent) in the CDPR scenario, and by the amount of 
rotation in the cooling flow models. 

In any case, the current explanations of the observational results are unsatisfactory In 
fact, Pellegrini et al. (1996) pose unsolved problems, especially concerning the explanation 
of the scatter in the Lx — Lb diagram, as a function of the flattening and rotation. It is 
found that Lx/Lb does not correlate with rotation and/or flattening. On the contrary, a 
very strong segregation effect is found: only galaxies with small flattening and low v rot /o~ c 
show a scatter in L X /L B , while flat objects and strongly rotating ones are systematically 
found only at low Lx/L-q. The main factor in determining the Lx/L-q segregation (and 
scatter, as a consequence) could be related to the slope of the galaxy surface brightness 
profile in the very center, as measured by HST observations (Pellegrini 1997, and references 
therein). This would tell us that a nuclear origin for the X-ray behavior of early-type 
galaxies should be taken into proper account (Ciotti & Ostriker 1997). It is remarkable 
how the independently chosen separation of early-type galaxies in cuspy and core galaxies, 
a separation operated using only photometric criteria, reproduces exactly the segregation 
in the Lx — Lb plane: while core galaxies show the entire scatter in Lx at fixed Lb, cuspy 
galaxies are invariably found only at low Lx- 

Other problems affect both the CDPR models and the cooling flow models. For 
example, the absence of massive cold disks in Es seems to support the CDPR scenario, 
while the claimed very low metallicity in the hot gas is against it, because a low metallicity 
contrasts with the assumption of a non negligible value of $sn. In fact, in principle, 
constraints on the SNIa rate can be given by estimates of the iron abundance in galactic 
flows (see, e.g., [Renzini et al. 1993| , and references therein). Under the assumption of solar 



abundance ratios, the analysis of the available data suggests a very low iron abundance, 
consistent with no SNIa's enrichment and even lower than that of the stellar component 
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(Ohashi et al. 1990; Awaki et al. 1991; Ikebe et al. 1992; Serlemitsos et al. 1993; 
Loewenstein et al. 1994; Awaki et al. 1994; Arimoto et al. 1997; Matsumoto et al. 1997). 
This is rather puzzling, in that a value as high as $sn = 0.3 agrees with the current 
optical estimates of the SNIa rates. However, some authors have found that more complex 
multi-temperature models with higher abundance give a better fit to the data (Kim et 
al. 1996; Buote & Fabian 1997). So, there are arguments on both sides of the abundance 
question, which is far from being a closed issue. In any case it is clear that the CDPR 
scenario is ruled out in case a null SNIa rate is clearly determined. 

An interesting feature of our models, that could be used to discriminate between the 
CDPR and the cooling flow scenario, is the spontaneous occurrence of the flow decoupling. 
This is at the origin of the striking feature of the transient cold filaments, which - if caught 
during their formation - should be observed first as soft X-ray emitting structures, and 
subsequently as Ha filaments. In model Al this takes place at tfu ~ 2 Gyr, while in model 
CI (the model with 6 = 4 kpc but similar to Al in all other input parameters) it occurs 
at tm ~ 10 Gyr. Due to this sensibility to the galaxy structure, it is plausible that for a 
few galaxies the decoupling of the flow takes place at the present epoch, and so these Ha 
features could be observed^. Indeed, it is tempting to identify the transient cold filaments 
with the soft X-ray clumps ( [Kim fc Fabbiano 1995| ) or Ha filaments ( |l'rinchieri fc di Serego 



Alighieri 1991| , Trinchieri, Noris, & di Serego Alighieri 1997) recently observed. We note 
how the size of the numerical and observed structures (of the order of ~ 10 kpc) are 
surprisingly similar. This could be an alternative explanation to the galactic drips proposed 
by Mathews (1997). 



4 Note that a flow decoupling - due to different physical reasons - is shown also in spherical 
power-law galaxies with low dark matter content and low SNIa heating ( Pellegrini fc Ciotti 



1997|) . 
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A final comment is on the X-ray surface brightness distributions of the computed 
models. As already found by BM for cooling flow models, rotation is the key factor in 
determining whether Ex is more or less round than E*, even in the presence of SNIa 
heating: rotation produces a flatter E x . A favourable feature of the CDPR scenario is 
that outflowing galaxies have rounder X-ray isophotes, which alleviates the problem of 
excessively flat Ex profiles, affecting our cooling flow models as well as the models of KM 
and BM. A detailed observational analysis of the relationship between optical and X-ray 
surface brightness for the E4 galaxy NGC 720 and SO galaxy NGC 1332 has been carried out 
by Buote & Canizares (1997, and references therein). While data for NGC 1332 are more 
uncertain, an analysis of the X-ray isophotes of NGC 720 shows that they are misaligned 
with and rounder than the optical ones. The authors conclude that the shape of E x is due 
to the flattening of the dark halo, which is higher than that of the stellar component. In 
fact, our simulations shows that, if rotation is the origin of flattening of the X-ray isophotes, 
then a decreasing flattening of Ex with galactocentric radius is expected, at variance with 
the observations. On the other hand, in our simulations flatter dark halos produce (during 
the inflow phase) flatter X-ray isophotes, as expected. 



5. Conclusions 

We have performed 2D numerical simulations of gas flows in SO galaxy models with 
different flattening, following the same recipe for the input physics adopted in CDPR, and 
exploring the effects of different amounts of ordered rotation and SNIa energy injection 
rates, and structural parameters. As a general rule, the WOI sequence which occurs in the 
ID models is replaced by a wind phase followed by a decoupled inflow. The main results 
can be summarized as follows: 

1. As expected models without SNIa heating (x = 4) host a cooling flow from the 
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beginning. The flow is considerably ordered, and essentially all the gas lost by the stars 
is cooled and accumulated in the galactic center. If rotation is present, the cold material 
settles in a disk on the galactic equatorial plane. 

2. Models with SNIa heating and x — 1 at the present epoch, invert their flow much 
earlier than the spherical ones with the same energy budget. In particular, as the SNIa 
rate decreases, the flow tends to revert its motion starting from the galactic polar regions. 
In this phase of strong decoupling, cold filaments are created at the interface between 
inflowing and outflowing gas. This strongly unstable phase lasts ~ 2 Gyr, after which the 
filaments are definitively accreted on the galactic center. Later on, the turbulences are 
unable to create regions of very cool, dense gas. The temperature may however fluctuate 
and gas cooler than X-ray temperatures is present. In the rotating model, the presence of 
SNIa greatly reduces the formation of the cold disk. In general, less cold matter is present 
in these models (~ 10 10 M ) than in those with x = 4 10 11 M ). 

3. Models with realistic values of the dynamical quantities are characterized by low 
X values. While for model Dl (b = 5 kpc, x — 0.7) the flow is always escaping without 
displaying any interesting feature, in model CI (b = 4 kpc, x — 0.85), the flow inversion 
occurs at t cc = 5.5 Gyr, and filaments form at tm ~ 10.5 Gyr. Of course, no cold mass is 
present in the first model; in the second, 5 10 9 M of cold gas form a small disk. 

4. For any fixed value of $sn, rotating models are X-ray fainter than the non-rotating 
ones, for the presence of the cold disk as well as for a different distribution of the hot gas. 
When compared to ID models, 2D models are much fainter in the X-rays. The former 
flows, in fact, have Lx similar to or even higher than Lsn after the flow inversion; the latter, 
instead, always have Lx < £sn- This is due to several combined factors: the loss of a SNIa 
energy fraction which is emitted by the gas below the X-ray temperature in high density 
regions (which are much larger than in spherical models); the loss of a SNIa energy fraction 
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used by the gas still outflowing from the external regions of the disk; the corresponding 
gravitational energy not released by this escaping material, and, for rotating models, the 
loss of thermalization of the gravitational energy of the cold disk. Another main difference 
of the temporal evolution of Lx is its noisy behaviour in 2D models. This is due to the 
presence of turbulence and to the formation of cold filaments. 

5. We have shown that the flattening of the galaxy in the CDPR scenario allows for 
a decoupling of the flow, in which different flow regimes are present in different part of the 
galaxy at the same time. We have found also that the flow decoupling is a general feature, 
i.e., all our models, except those without SNIa heating and the model with x substantially 
lower than unity, show a flow decoupling. We note however that the Lx associated to a 
decoupled flow when x < 1? is n °t high: it is true, on the contrary, that the flows are 
decoupled even when % > 1, as the galaxy loses gas from the disk even with a low SNIa rate. 
So, the energetical analysis carried out by CP is shown to be qualitatively verified: objects 
with a low x value are found preferentially at low Lx, and at variance with the cooling 
flow models the effect of the flattening is shown to be important in determining the flow 
evolution. 

6. For non- rotating, inflow models, the edge-on Ex is rounder than the corresponding 
E* of the parent galaxy; when rotation is present, Ex becomes flatter, and very similar to 
E*. For models with SNIa heating, the E x morphology evolves with time, from a smooth, 
peanut-shaped geometry to a generally rounder shape than the stellar profile, through 
disturbed phases during the presence of cold filaments. When rotation is present, in the 
galagtic center the Ex isophotes are somewhat flatter than E*, at variance with non-rotating 
models in which Ex remains more spherical up to the center. As a general rule, even for 
low line-of-sight inclination angles over the galactic equatorial plane, the isophotes of E x 
and E* tend to become similar. 
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All these new features considerably enrich the old CDPR scenario, and represent a step 
further in the direction of a more realistic simulation of a situation that is proving to be 
very complex. In future, we plan to perform numerical simulations of gas flows in models 
tailored on observed X-ray emitting SO galaxies, in order to improve our understanding of 
the much richer phenomenology of aspherical gas flows. 
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Fig. 1. — The gas density and velocity fields of model Al at four different times (0.46, 2.16, 
3.37, 15) Gyr (top left, top right, bottom left, bottom right). In the top left panel the dashed 
lines represent the stellar density distribution. In the second panel the large high density 
transient structure is apparent. Note how the velocity field in the second and third panels 
is dominated by large vortices, while in the fourth panel the velocity is very low, and nearly 
random. 

Fig. 2. — The gas temperature and velocity for model Al (high resolution simulation) at 
t ~ 2 Gyr, when the cold filament is completely developed. Note how the disk is still strongly 
degassing. 

Fig. 3. — The time evolution of the mass and luminosity of gas flows for the models discussed 
in the text. Panel (a) illustrates the behaviour of the hot gas of model Bl (t?g N = 0, k — 1; 
dashed line), and of model B2 ($ S n = 0, k — 0; dotted line). The solid line represents the 
total mass lost by the stars. Panel (c) represents the L x evolution for the same models. 
Panel (b) shows the mass budget for the models with -# SN = 0.3. Dotted lines represent the 
evolution of the hot gas mass for model Al (k = 1, line 1); A2 (k = 0, line 2); CI (k = 1, 
line 3); Dl (k — 1, line 4). Dashed lines show the total gas mass inside the galaxy. The solid 
line again represents the total mass lost by the stars. Panel (d) shows the evolution of Lx 
for these models (dotted lines), and the solid line represents Lg N . 

Fig. 4. — The optical (solid lines) and X-ray (dotted lines) surface brightness distributions 
for model Al, at the same representative times as in Fig. 1 (time increases from top to 
bottom). The optical isophotes are separated by one magnitude. The X-ray isophotes are 
logarithmically spaced: the lowest value (-8.51 in cgs units) is the same in all panels, the 
maximum is -3.34 (a), -2.97 (b), -3.37 (c), -4.47 (d). 

Fig. 5. — The gas temperature and velocity for model CI at t ~ 10 Gyr, when the cold 
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filament is completely developed. 
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Table 1. Models parameters. 
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a All models have L B = 10 10 L Q , M* = 10 11 M Q , and M h /M» = 2 
b in kpc 

c three-dimensional central velocity dispersion in km s _1 
d maximum rotational velocity in km s _1 

°for the rotating models computed under the assumption of 
complete thermalization of ordered motions 
f in Gyr 
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